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Abstract. We extend the analysis of the renormalization group flow in the two-dimensional Hubbard 
model close to half-filling using the recently developed temperature flow formalism. We investigate the 
interplay of d-density wave and Fermi surface deformation tendencies with those towards d-wave pairing 
and antiferromagnetism. For a ratio of next nearest to nearest neighbor hoppings, t'/t — —0.25, and 
band fillings where the Fermi surface is inside the Umklapp surface, only the d-pairing susceptibility 
diverges at low temperatures. When the Fermi surface intersects the Umklapp surface close to the saddle 
points, d-wave pairing, d-density wave, antiferromagnetic and, to a weaker extent, d-wave Fermi surface 
deformation susceptibilities grow together when the interactions flow to strong coupling. We interpret 
these findings as indications for a non-trivial strongly coupled phase with short-ranged superconducting 
and antiferromagnetic correlations, in close analogy with the spin liquid ground state in the well-understood 
two-leg Hubbard ladder. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 74.72.-h High-Tc compounds 
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1 Introduction 



The two-dimensional (2D) Hubbard model is a focus of 
electron correlation theory because of its relevance to high- 

■ temperature superconductivity [EJ and other correlation 
phenomena such as metal-insulator transitions and itin- 

, erant ferromagnetism H . Apart fromproposed exotic spin- 
liquid or non-Fermi liquid states a widely adopted 
approach is to classify possible ground states of the model 

. in terms of broken symmetries. This can be done for the 

' strong-coupling versions of the Hubbard model by invok- 
ing spin-charge separation |J or large-spin approaches . 
But for moderate interaction strengths it is more com- 
mon to invoke a mean field description of long-range or- 
dered ground states with anomalous expectation values 
which are bilinear in electron creation and annihilation 
operators. These theories are either constructed with the 
bare Hubbard interaction ®^ or with an effective inter- 
action which is postulated or calculated perturbatively, 
e.g. by one-loop renormalization group (RG) treatments 
or parquet summations. Typically these treatments are re- 
stricted to the weakly correlated case, but often they can 
give qualitatively satisfactory results such as the antiferro- 
magnetic state at half-filling or a <i a; 2_ J( 2-superconducting 
ground state at a reduced band filling in the 2D Hubbard 
model 
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Recently two independent proposals enlarging the vari- 
ety of potential symmetry-broken states have been made. 
One possibility was raised by Halboth and Metzner jnj, 
based on an analysis of the Landau function extracted 
from a RG flow. They pointed out that the system might 
become unstable against a Pomeranchuk instability, i.e. 
a spontaneous deformation of the Fermi Surface (FS) re- 
ducing its symmetry from square to orthorhombic. Qual- 
itatively similar results were obtained by Valenzuela and 
Vozmediano fl~5|| in a mean- field study of an extended Hub- 
bard model. However with this technique and also with the 
Wilsonian RG formalism based on an infrared (IR) cutoff 
used in Ref . 1 1 it is difficult to compare the strength of 
the FS deformation tendencies with other potential singu- 
larities, e.g. in the AF or d-wave superconducting channel. 
Here we use the recently developed temperature-flow for- 
malism [ fl6| , to analyze these FS deformation tendencies 
together with other channels in an unbiased way. 

A different t ype of symmetry breaking was consid- 
ered by Nayak fll7f and Chakravarty et al. |Q. They 
proposed that a d-density-wave state (closely related to 
flux phases jl9| or orbital antiferromagnets pQ]) forms at 
a high temperature in the underdoped high-T c cuprates 
and is responsible for the many anomalous properties like 
the notorious pseudogap. This proposal has been contro- 
versial for a number of reasons pi], p2|, not least because 
of the absence of experimental evidence that the onset of 
the pseudogap is accompanied by an actual phase tran- 
sition rather than a continuous crossover. The onset of a 
phase with broken translational and time reversal symme- 
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tries should cause anomalies in thermodynamic and other 
properties which have not been observed. Nonetheless it 
is interesting to investigate whether d-density-wave ten- 
dencies play a significant role in the RG flow to strong 
coupling. 

In this paper we show that in the one-loop RG tem- 
perature-flow to strong coupling in the t-t' Hubbard model 
at U ~ 3t both the FS deformation and d-density wave 
tendencies grow. However their growth is weaker than that 
of d-wave pairing and AF fluctuations, so that they do 
not appear as the dominant instabilities of the model. 
Nonetheless, similar to our findings with a RG IR-flow 
in Ref. 
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we find a particular density region - the so- 
called saddle-point regime which occurs when the FS is 
close to the saddle points of the band dispersion at (n, 0) 
and (0, 7r) - where d-wave pairing, AF, d-density wave and 
d-wave Pomeranchuk fluctuations are intimately related 
and grow together in the range where the one-loop flow 
is credible. This contrasts with the behavior for fillings 
smaller than the van Hove filling where only the d-wave 
pairing correlations continue to grow towards lower tem- 
peratures, while the others are cut off below some tem- 
perature. 

The qualitative picture obtained here is a confirma- 
tion and extension of the results described in Ref. fll2|] . 
The main motivation for the present paper is twofold: a) 
we show that the new temperature-flow RG scheme re- 
produces the results of the IR-flow RG Jl^] and yields an 
even clearer picture; b) we now also include d-wave FS de- 
formations and d-density wave tendencies which were not 
considered in the previous analysis. 

In the following we first the describe the calculational 
scheme and the symmetry-breaking tendencies which are 
investigated. Then we discuss our numerical results for the 
case of zero and non-zero values of the next-nearest hop- 
ping t' and conclude with a comparison of our observations 
to ladder systems, which we argue give some insights into 
the nature of the strongly coupled state. 



2 The calculational scheme 



The Hamiltonian for the 2D t-t' Hubbard model is 



U 



with onsite repulsion U and hopping amplitudes t and t' 
between nearest neighbors (n.n.) and next-nearest neigh- 
bors (n.n.n.) on the 2D square lattice. We apply the so- 
called temperature-flow RG scheme introduced recently Jl6| 
in a A-patch implementation that covers the full Fermi 
surface . Similar to the approaches in Refs. [I4 P|| 

the T-flow scheme is derived from an exact RG equation. 
However, a low energy IR cutoff is not used, and instead 
the temperature T itself is used as the flow parameter. 
This new formulation allows an unbiased comparison be- 
tween AF and FM tendencies, whereas RG schemes with 
a flowing IR cutoff artificially suppress particlc-hole exci- 
tations with small wave vectors, e.g. long wavelength den- 
sity fluctuations. For a detailed discussion of this matter 



and a derivation of the T-flow RG equations the reader is 
referred to Ref. [[l6) . 

The RG formalism yields a hierarchy of differential 
equations for the one-particle irreducible n-point vertex 
functions, Tj" , as functions of the temperature. Integra- 
tion of this system of equations gives the temperature- 
flow. As initial condition we assume that at a high tem- 
perature To, the single-particle Green's function of the sys- 
tem is simply Go(iuj, k) = [iuj— e(k)] -1 and the interaction 
vertex is given by a local repulsion, U . This is justified if 
To is sufficiently large, as wavevector-dependent perturba- 
tive corrections to selfenergy and four-point vertex decay 
at least oc 1 / Ton. We truncate the infinite system of equa- 

. . (n) 

tions by dropping all vertex functions T^, with n > 4. 
In the present treatment we also neglect selfenergy cor- 
rections and the frequency dependence of the vertex func- 
tions. This restricts the scheme to one-loop equations for 
the spin-rotation invariant four-point vertex T^ . Starting 
with weak to moderate interactions, we follow the T-flow 
of as T decreases. 



The four-point vertex is deter mined b; 



a cou- 
pling function V T (k 1: k 2 , k 3 ) (see Refs. The 
numerical implementation follows the work of Zanchi and 
Schulz jlo| and was already explained in Refs. jl2] and 
[ pi] . We define elongated phase space patches around lines 
leading from the origin of the BZ to the (±7r, ±7r)-points, 
and approximate Vr(ki, k.2, k 3 ) by a constant for all wave 
vectors in the same patch. We calculate the RG flow for 
the discrete subset of interaction vertices with each FS 
patch represented by a single wave vector and with the 
initial condition Vt (ki, k2, k 3 ) = U. Most calculations 
were performed using 48 patches. Furthermore we calcu- 
late the flow of several static susceptibilities, as described 
below. In addition we consider the flow of the couplings 
h c (k) of quasiparticlcs at different positions on the Fermi 
surface to uniform static charge fields. This allows us to 
analyze which classes of coupling functions and which sus- 
ceptibilities become important at low T. For a large pa- 
rameter range, we observe a flow to strong coupling, i.e. at 
sufficiently low temperature some components of the cou- 
pling function Vr(ki, k2, k 3 ) become larger than the band- 
width. The approximations mentioned above fail when the 
couplings are too large. Therefore we stop the flow when 
the largest coupling exceeds a high value larger than the 
bandwidth, e.g. Vr,max — 18i. This defines a characteristic 
temperature T* of the flow to strong coupling. 



3 Antiferromagnetic, c?-wave pairing, 
d-density waves and d-wave Pomeranchuk 
fluctuations 

In this section we introduce fermionic coupling terms to 
order parameters or static external fields corresponding 



1 A wavevector-independent contribution can be absorbed 
into the chemical potential fi. Flows with fixed particle number 
instead of fixed fi give qualitatively similar results. 
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to different symmetry breaking channels. Within the RG 
treatment described in the previous section, we then calcu- 
late the one-loop renormalizations of these coupling terms 
and thus obtain information on the growth of the cor- 
responding fluctuations in the low-temperature state. A 
divergence of one of these external couplings signals a 
crossover to a region with a strong tendency to ordering 
which could become an actual finite temperature phase 
transition when coupling in a third spatial direction is 
added. 

First we introduce the ci-wave pairing field <? S! d_ sc = 
— <£_ SiC j-sc, which couples to the electrons via 



(p 



s,d—sc 



k 



=(k)4. 



(1) 



with a coupling constant hd- sc Qc). As initial condition 
at T = To we assume a d x 2_ y 2 form factor hd— S c(k) = 
(cosfcz — cos k y ) /\/2. Like the other coupling constants 
defined below, hd~ S c will develop a k-dependence through 
the perturbative corrections in the course of the temper- 
ature flow towards lower T and higher harmonics of the 
same representation of the point group will be admixed. 
Next we also define a (s-wave) spin-density wave field 
«&AF- Due to spin-rotational invariance it is sufficient to 
consider staggered fields in the spin quantization direc- 
tion, 



<? z ,af ^/iAF(k) (cl +Q s c K , 



k+Q 



-.,Ck,-; 



(2) 



Here the momentum transfer Q = (tt, n) between created 
and annihilated particles corresponds to an alternating 
field acting on the electron spins on n.n. sites. The initial 
condition at high temperatures is /iaf (k) = 1. The pro- 
posed Pomeranchuk deformations also have d x 2_ y 2 sym- 
metry, the coupling term is given by 



hm d _ P \ hd- 



k+q.s^k.s ' 



(3) 



with initial condition /id-p(k) = (cosfc^ — cosk y ) /\/2 at 
T . A non-zero d-w&ve Pomeranchuk field ^d-p lifts the 
degeneracy of the saddle points at (ir,0) and (0,7r) and 
breaks the fourfold symmetry of the kinetic energy. The 
(i-density wave fluctuations couple to particle-hole pairs 
with momentum transfer Q, 



d-DW 



frd-dw(k) 4 



k+Q,s 



(4) 



with initial conditions chosen as 



M-DW 



(k) = (cos k x — cos k v 



at the initial temperature To. In the one- loop T-flow RG 
treatment 1 16, 2j| the coupling constants h are then renor- 
malized according to 



dT 



(k) 



1 



(k')y T (k,-k,k') 



^AF(k) = ~J2 hAF ^ VT ^ k '' k ' + Q ) 
' ' L k' 



dT 



Tp H (k',k' + Q), 



(5) 



(6) 



dT 



M-DWl 



dT 



k ) = -JT E^-Dw(k')^ip H (k',k' - Q) 

k' 

• [V T (k, k',k' + Q)- 2Vr(k, k', k + Q)] , (7) 
W = E ^-P( k ') ^^(k', k') 

k' 

■[y T (k,k',k')-2y T (k,k / ,k)] . 



(8) 



The expressions on the right hand side are evaluated in ab- 
sence of the external fields. Nl denotes the number of lat- 
tice sites and the particle-particle and particle-hole loops 
are given by 



i P p(k,k') = T]r — 



i 



iuj n - e k -iw n - e k / 



and 



L PH (k,k') = rV — 



iu> n - e k iu n - e k < 



where ioj n denotes the fcrmionic Matsubara frequencies 
u> n — irT(2n + 1) with n = 1,2,.. .. The susceptibilities 
for the various channels are also renormalized by temper- 
ature derivatives of bubble diagrams with the correspond- 
ing couplings at the vertices. For example for the d-wave 
pairing susceptibility we obtain 



dT 



Xd- 



— ^fr d „ sc (k) 



dT 



L PP (k,-k)^_ sc (-k). (9) 



Analogous expressions hold for the other susceptibilities. 
As initial condition we assume that all susceptibilities are 
zero at the initial high temperature, T . 



4 Results for the 2D Hubbard model 

Here we describe our numerical results obtained with the 
A-patch implementation of the one-loop T-flow RG scheme 
with N — 48 patches. For all cases discussed below, the 
RG flow goes to strong coupling. This means that some 
components of the coupling function become larger than 
the perturbative range when the temperature is lowered 
far enough. We define a characteristic temperature T* for 
the flow to strong coupling determined as the temperature 
when the largest coupling reaches V max = 18*. 
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4.1 The case t' = 

First let us discuss results for vanishing next-nearest neigh- 
bor hopping, t' — 0. The flow of the susceptibilities is 
summarized in Fig. |l| 

For the square Fermi surface at half band filling, i.e. 
fj, = 0, we find that the AF susceptibility is the lead- 
ing divergent susceptibility, consistent with the expected 
AF ordered ground state. We also find weaker growing d- 
wave pairing and d-density wave susceptibilities. At = 0, 
these two channels have exactly the same temperature 
dependence due to the particle-hole symmetry; e(k) = 
— e(k + Q). For finite doping, i.e. fi ^ 0, the flow splits this 
degeneracy and the d-wave pairing generally grows more 
strongly than the d-density wave channel. As the electron 
density is further reduced below half-filling, the nesting in 
the (tt, 7t) particle-hole channel is increasingly weakened. 
Correspondingly, at a critical doping that depends on the 
interaction strength, the growth in the AF and d-density 
wave channels get cut off at low T and the d-wave pairing 
is the only divergent channel that remains. This crossover 
from the AF dominated flow to the d-wave pairing domi- 
nated instability for t! = is very similar to the results of 
Zanchi and Schulz |l0| and Halboth and Metzner jll]] who 
used a TV-patch scheme with flowing IR cutoff. Within the 
latter schemes it is difficult to compare the flow of small-q 
particle-hole susceptibilities, e.g. Pomeranchuk FS defor- 
mations, with the flow of AF or superconducting suscepti- 
bilities. There is no such difficulty in the temperature flow 
scheme. Thus we can directly compare the Pomeranchuk 
channel with the above-mentioned other tendencies. The 
main conclusion arising from the data shown in Fig. |l| is 
that the d-wave Pomeranchuk tendencies grow somewhat 
to lower temperatures, but apparently they do not repre- 
sent the leading instability for the t' = Hubbard model 
with only a local onsite repulsion as in initial interaction. 
In this paper we do not give a systematic analysis how the 
leading instability may change when the initial interaction 
is varied, e.g. by including a n.n.n. repulsion V or a spin 
exchange coupling J. This question was addressed by Binz 
et al. [£6| f° r t ne half-filled case and U — > 0. These authors 
found a d-density wave regime for sufficiently large posi- 
tive n.n. Heisenberg coupling J between an AF regime at 
small n.n. repulsion V and s-charge density wave regime 
at AV > U. We find a similar crossover between the two 
latter regimes for U > 0.3i. However a RG flow where the 
d-density wave susceptibility grows most strongly was only 
found for a limited parameter region around J > 0.7C/ and 
AV « 0.7U. 



4.2 Moderate n.n.n. hopping t' = -0.25t 

Next we focus on the case of n.n.n. neighbor hopping 
t' = — 0.25t and band fillings n < 1 for which the FS 
passes close to the saddle points. Upon increasing the par- 
ticle density through the van Hove filling at which ijl lies 
exactly at the van Hove singularity, i.e. /i = At', the FS 
starts to intersect the Umklapp surface which connects 
the (tt, 0) and (0, 7r) points with straight lines. As we will 




* ~" i o — * — ><~ ~ 

0.5 1 0.5 1 

T/t T/t 



Fig. 1. Temperature flow of AF (thick solid line), d- wave sc 
pairing (thin solid line), d-density wave (dashed-dotted line) 
and d-Pomeranchuk susceptibilities for different values of the 
chemical potential and t' = 0. The cross and the triangle at 
the T axis denote the temperatures where the largest couplings 
exceed 5i and Wt. Half filling corresponds to /i = 0. 



show, this gives rise to a remarkable qualitative change 
in the flow. Two typical shapes of the the non-interacting 
Fermi surfaces for densities below and above the van Hove 
filling arc shown in Fig. |^. 

In Fig. ^ we show the evolution of d-wave pairing, AF, 
d-density wave and d-wave Pomeranchuk susceptibilities 
for different values of the chemical potential as the tem- 
perature is lowered. In the first upper two plots for which 
[i < At' , we observe that only the d-wave pairing suscepti- 
bility grows rapidly at low temperatures T — + T* > 0. The 
other susceptibilities reach a maximum at a T max > T* 
and decrease again as T — »• T*. For the next two plots 
the FS intersects the Umklapp surface close to the saddle 
points and the flow is qualitatively different: now all four 
susceptibilities grow as T — > T* . As described in Ref . |Q , 
the flows of AF and d-wave pairing susceptibilities are very 
similar for a broad density region above the saddle point 
filling and it appears to be implausible that the strong 
coupling state is a simple symmetry-broken state corre- 
sponding to a single type of ordering. More than that, 
the d-density wave and d-wave Pomeranchuk susceptibil- 
ities also grow as T — > T* . But for our parameters these 
two channels do not constitute the leading instabilities. 
Nevertheless they seem to be part of a common mecha- 
nism arising from the special location of the FS close to 
saddle points which causes mutual reinforcement of dif- 
ferent channels, most prominently AF and d-wave pairing 
described in Ref. [l^ . The main reason for this effect is 
that the scattering processes connecting the two inequiva- 
lent saddle point regions involve particle pairs with small 
total momentum that scatter with wavevector transfers 
close to (tt, n). Therefore these scattering processes are 
strongly driven by both the Cooper and the (tt, tt) particle- 
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Fig. 2. Non-interacting Fermi surfaces with t' = — 0.25t for 
/i = — 1.05t (left plot, (n) ~ 0.75 per site) and /i = — 0.95t 
(right plot, (n) f» 0.82 per site). The dots denote the location 
of the 48 points for which the flow of the coupling function is 
evaluated. The dashed square denotes the so-called Umklapp 
surface. 




0.5 1 0.5 1 0.5 1 
T/t T/t T/t 



Fig. 3. Temperature flow of AF (thin solid line), d-wave pair- 
ing (thick solid line), d-density wave (dashed-dotted line) and 
d-Pomeranchuk susceptibilities for different values of the chem- 
ical potential and t' = — 0.25i. The cross and the triangle at the 
T axis denote the temperatures where the largest couplings ex- 
ceed 5t and lOt. Band fillings less than the saddle point filling 
correspond to n < —t. (n) « 0.93 per site for fx — — 0.75t. 



hole channel and this causes a mutual reinforcement of the 
fluctuations generated by these channels. Thus it is inter- 
esting to speculate that the strong coupling state will em- 
body all these channels as short range correlations, instead 
of selecting a single ordering channel and suppressing the 
others. Indeed, as discussed in the concluding section, ex- 
amples for such ground states are known from ladder sys- 
tems. Furthermore, other numerical investigations of the 
doped 2D t-J model showed the existence of d-density 
wave correlations in Gutzwillcr-projected d-wave pairing 
variational wave- functions |24| and in the lowest energy 
d-wave paired states in exact diagonalization on a 32-site 
cluster {§. 

We note that if we continue to follow the one-loop 
RG flow for the 2D system which ignores selfenergy cor- 
rections, the common growth of the different channels 




0.1 1 2 0.1 1 2 0.1 1 2 

T/t T/t T/t 



Fig. 4. Temperature flow of AF (thin solid line), d-wave pair- 
ing (thick solid line), d-density wave (dashed-dotted line) and 
d-Pomeranchuk susceptibilities for different values of the ini- 
tial interaction U and /i = — 0.92t at t' = — 0.25i. The cross 
and the triangle at the T axis denote again the temperatures 
where the largest couplings exceed 5t and lOt. 

does not continue to arbitrarily low temperatures. Typi- 
cally a decoupling takes place once the overlap between 
the channels becomes too small due to the restriction 
of the relevant scattering processes at low temperatures 
to thinner and thinner neighborhoods of the Fermi sur- 
face. Signs for this are visible for example in the flow of 
d-density wave susceptibility in the last plot of Fig. 
However we are convinced that one should concentrate on 
the flow down to temperatures where the growing inter- 
actions reach values comparable to the bandwidth. This 
is also the scale where the corrections to the one-loop 
RG equations such as one- and two-loop selfenergy ef- 
fects start to become non-neglig iblc @,|l§]. In the saddle 
point regime described above the decoupling of the differ- 
ent channels only occurs at lower temperatures where the 
one-loop flow without selfenergy corrections has lost its 
validity as the couplings are much larger than the band- 
width. Obviously, the density range where this common 
flow of several channels occurs becomes larger with in- 
creasing interaction strengths U. The reason is that at 
higher temperature scales a wider shell around the FS 
participates in the RG flow and the overlaps between the 
scattering processes driving the individual fluctuations are 
larger. This is illustrated in Fig. ^, where we plot the flow 
of the susceptibilities for the same chemical potential but 
for different initial interactions U. While for U = 2.5t the 
flows of AF, d-density wave and d-wave Pomeranchuk sus- 
ceptibilities go through a maximum before the couplings 
exceed the bandwidth, for U > 3i their flow grows mono- 
tonically over the range where the one-loop flow applies. 
Thus we expect that there is a nonzero minimal interac- 
tion strength above which the coupled saddle point flow 
determines the ground state properties. For smaller inter- 
actions the low energy physics may be dominated by a 
single channel. 

With the new T-flow RG scheme we again find an 
anisotropic suppression of the charge compressibility in 
the saddle point regime, similar to our findings in Ref. p2| . 
This is illustrated in Fig. [jj, where we plot the T-flow of the 
couplings to external static uniform charge fields, h c (k). 
These quantities are calculated in a way analogous to the 
Pomeranchuk coupling constants hd-p(\i), with the differ- 
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Fig. 5. Temperature flow of couplings /i c (ki?) of quasiparti- 
cles on different points kp on the FS to uniform static charge 
fields. The insets show the flow normalized to the RPA values 
for the corresponding temperatures. The solid (dashed) lines 
correspond to k_F points close to the saddle points (BZ diago- 
nals). The left plot is for /i = ~l.lt (band filling less than the 
van Hove value, d-wave regime, band filling less than van Hove 
filling), and the right plot for /i = — 0.95t (saddle point regime, 
band filling slightly larger than van Hove filling). 



ence that the d-wave form factor in Eq. (||) is replaced by 
unity. In the left plot the density is less than the van Hove 
filling corresponding to the pure d-wave regime and the 
charge couplings h c (k) are suppressed uniformly around 
the FS. Moreover, as shown in the inset, the RG flow of the 
interactions enhances the charge couplings h c (k) with re- 
spect to isotropic RPA (random phase approximation) val- 
ues h° c calculated without RG flow of the interactions. In 
the right plot, we show the same quantities for \i = — 0.95i, 
corresponding to the saddle point regime at band fillings 
slightly larger than the van Hove value. Here the suppres- 
sion of the charge couplings is strongest in the FS region 
close to the saddle points. Furthermore, compared to the 
RPA values, the charge couplings for this FS region are 
strongly decreased at low T by the RG flow of the inter- 
actions, while they increase near to the BZ diagonals. We 
emphasize that these tendencies are rather weak because 
we do not let the interactions become too strong. Never- 
theless, our findings with the new T-flow RG scheme are 
again consistent with the hypothesis of a Fermi surface 
truncation [^,|l^l at the saddle points where the scatter- 
ing processes grow most strongly, with remaining FS arcs 
around the BZ diagonals. We believe that the effects ob- 
served as tendencies in the weak-coupling approach will 
be strongly amplified in a full strong coupling description. 
We add that a similar suppression of the k-space resolved 
compressibility at the saddle points is observed within a 
Landau-Fermi liquid framework p8f and was also found 
in quantum Monte Carlo calculations for the half filled 
Hubbard model by Otsuka et al. f p9| . 

Lastly we comment on the behavior of the T-flow RG 
scheme as the density is increased towards half filling. The 
saddle points move to larger negative band energies and 
the coupling between the channels becomes less effective 
(see middle and right lower plots in Fig. ||). The AF sus- 
ceptibility grows most strongly at higher T. At lower tem- 
peratures the T-flow RG scheme exhibits a stronger sen- 
sitivity to the imperfect nesting introduced by the n.n.n. 
hopping term t' than in the IR-flow RG scheme we used 



previously pig . This shows up in a saturation of the AF 
susceptibility flow for weaker initial interactions at low T 
before the couplings have reached values larger than the 
bandwidth. At half- filling an initial value of U > 3.5£ is re- 
quired to obtain the dominant AF phase found in Ref. [fi~2| . 



5 Conclusions and comparison with the 
two-leg Hubbard ladder 

We have investigated the flow of d-wave pairing, spin den- 
sity wave, d-density wave and d-wave Pomeranchuk sus- 
ceptibilities in the two-dimensional Hubbard model near 
the van Hove filling. For a n.n.n. hopping t' = — 0.25i we 
found two distinct regimes. For densities well away from 
the van Hove density, only a single susceptibility diverges 
at low T. For a range of band fillings n < 1 around, or 
somewhat above, the van Hove filling, whose width de- 
pends on the strength of the initial interaction, several 
susceptibilities grow together when the interactions flow 
to strong coupling. In particular, from the interplay be- 
tween d-wave pairing and AF tendencies near the saddle 
points, also described in Ref. ]12| ], we conclude that the 
strong coupling phase, in the saddle point regions, is not 
given by a dominance of a single ordering such as super- 
conducting, magnetic or charge ordering. Furthermore we 
again find tendencies towards incompressibility of the FS 
Pomeranchuk regions near the saddle points which can be 
interpreted as indication for a FS truncation. Although 
the d-density wave and d-wave FS deformation modes 
get enhanced at lower temperatures, but do not repre- 
sent the leading instabilities for the t-t 1 Hubbard model 
with weak to moderate repulsive onsitc interaction. This 
renders a spontaneous symmetry breaking in these chan- 
nels unlikely. Nevertheless there may be significant short 
ranged correlations of the respective types. 

A similar RG flow to strong coupling is also found 
in the half-filled two-leg Hubbard ladder. This quasi-one- 
dimensional system has been extensively analyzed both 
from the weak-coupling perspective |30|, pl| , ^2l a nd by nu- 
merical methods for stronger onsite repulsion []33[ . Through 
these studies it became clear that the ground state for all 
U is given by an insulating spin liquid, i.e. spin and charge 
fluctuations are gapful and decay exponentially. The FS 
is fully truncated although no symmetry breaking occurs. 

The weak-coupling analysis 0,^f) typically starts with 
a one-loop RG scheme which selects the relevant cou- 
plings. Subsequently the effective Hamiltonian, composed 
of the kinetic energy and the relevant couplings, is bosonized 
and the properties of spin and charge excitations can be 
deduced. Our comparison with the 2D Hubbard model 
focuses on only this one-loop RG step. In the one-loop 
g-ology for the half-filled two- leg ladder, one finds seven 
relevant coupling constants gi, i = 1, . . . , 7 which describe 
scattering processes between the 4 Fermi points, 2 in the 
even parity band and 2 in the odd band. In the RG flow 
with decreasing infrared cut-off A and initial condition 
9i(Ao) = U these couplings constants diverge asymptoti- 
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cally together as 



9M) 



log(A/A c ) 



(10) 



when A is reduced towards a finite critical scale A c . The 5? 
reach fixed ratios in the weak interaction limit. Due to the 
divergence of the coupling functions, the one-loop renor- 
malizations drive several couplings to external fields to in- 
finity. Like in the 2D system several types of fluctuations 
are driven by common scattering processes and mutually 
reinforce each other. As mentioned in the last section, in 
the 2D model we need a finite initial interaction strength 
U to observe the coupled flow of several channels for a 
reasonably wide parameter region. In contrast with that 
in the two-leg ladder the mutual reinforcement is present 
at all energy scales and occurs for arbitrarily weak inter- 
actions. More precisely, with the asymptotic flow given by 
Eq. (jic]) we find that in the half-filled ladder AF, d-wave 
pairing, d-density wave and spin density wave couplings 
grow with the same strengthn. Equivalently, if we write 
down the effective Hamiltonian of the system close to the 
critical scale, the effective coupling constants for mean- 
field decouplings in d-wave pairing, AF and d-density wave 
channel reach the same absolute value in the asymptotic 
flow very close to the instability. Since the effective den- 
sities of states for the three types of mean-field decou- 
plings are the same, their potential ground state energies 
are degenerate as well. Nonetheless, as we know from the 
bosonization |31| and numerical studies | p3| , the ground 
state does not exhibit (quasi-)long-range order of any type 
and the superconducting and magnetic correlations which 
seem to diverge in the one-loop treatment remain strictly 
short-ranged. 

As another similarity between the 2D system and the 
two-leg ladder we mention that also in the ladder the 
forward scattering amplitudes calculated within the RG 
scheme exhibit a d- wave- type structure with different signs 
for the forward scattering onto the same band and onto 
the other band. However a shift of the Fermi points in- 
duced by the flow of the forward scattering processes was 
found in a weak-coupling RG treatment only for the case 
of unequal Fermi velocities in the two bands ]34[ ]. 

Finally we add that in the 2D Hubbard model with 
t' < for band fillings n > 1 corresponding to electron- 
doped case the situation is different from the hole-doped 
case described above|p|. For n > 1 the (77, 0) and (0, tt) 
regions of the BZ, that are crucial for the channel coupling 
leading to the combined flow in the saddle point regime 
similar to the two-leg ladder flow, are at larger negative 
band energies and exert only little influence on the low 
temperature physics. Correspondingly, upon moving away 
from half filling by increasing n > 1 and suitable parame- 
ters, the flow to strong coupling rapidly changes from an 
AF regime with high T* to a <i-wave pairing regime with 
low T* without the occurrence of an intermediate saddle 
point regime. 



2 This is also true for the 2D model restricted to the saddle 
point regions, see Ref. [M 



In summary, the results of the recently developed T- 
flow one loop RG scheme applied to the t-t'-U Hubbard 
model in 2D, are qualitatively consistent with those we 
found earlier using an IR-flow RG scheme. The T-flow RG 
scheme has the advantage that the competition with long 
wavelength instabilities such as a d-wave FS Pomcranchuk 
distortion can be analyzed in an unbiased way. Our results 
show that this susceptibility and also the d-density wave 
susceptibility are never the dominant ones in a Hubbard 
with just an onsite repulsion, U although for densities in 
the saddle point regime they appear as divergent suscep- 
tibilities together with those towards d-wave pairing and 
AF order. These results agree with the proposal we made 
earlier jL2| regarding the nature of the strong coupling 
phase that appears at temperature T < T* . Based on the 
close analogy with the one loop RG results for the two 
leg ladder at half-filling, we proposed that at T w T* 
a crossover occurs to a strong coupling phase with a FS 
which is partially truncated in the vicinity of the saddle 
points by the formation of an insulating spin liquid which 
in turn is not characterized by any, even quasi, long range 
order. The remaining FS arcs centered on the BZ diago- 
nals are coupled in the Cooper channel to the insulating 
spin liquid and this coupling is the origin of a transition to 
d-wave superconductivity at a lower T < T* . Note we find 
that T* is a strong function of the density and increases 
from zero with increasing density as the van Hove filling is 
crossed. This behavior of T* is, we argue, consistent with 
the critical density and rising temperature scale found by 
Tallon and Loram J3f| . 

Of course our one-loop approach does not yield a con- 
trolled theory of the strongly coupled phase below the 
characteristic temperature T* and the second step in the 
weak-coupling analysis of the exemplary ladder systems, 
the bosonization of the effective Hamiltonian, cannot be 
performed in the 2D case. Thus we cannot prove the prop- 
erties of the 2D strong coupling phase for the time being. 
More powerful methods need to be developed in order to 
establish the nature of the strong coupling phase at tem- 
peratures T < T*, to test our proposal that it is closely 
analogous to that of the 2-leg ladder and also to examine 
possible connections to the gauge theories which proceed 
directly from strong coupling and also predict strong fluc- 
tuations coexisting in several channels such as the d-wave 
pairing and d-density wave channels p2j . 

We thank M. Sigrist, M. Troyer, D. Poilblanc, W. Met- 
zner, and C. Nayak for helpful discussions. C.H. acknowl- 
edges support by the Swiss National Science Foundation 
(SNF) and the Emmy-Noether program of the Deutsche 
Forschungsgemeinschaft (DFG). The numerical calcula- 
tions were performed on the Asgard Beowulf cluster of 
ETH Zurich. 
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